Carrier heating and high-order harmonics generation 
in doped graphene by a strong ac electric field 
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The nonlinear response of electrons (holes) in doped graphene on ac pumping is considered the- 
oretically for the frequency region above the energy relaxation rate but below the momentum and 
carrier-carrier scattering rates. Temporally-dependent heating of electrons by a strong ac field, 
which is described within the energy balance approach, leads to an effective generation of high- 
order harmonics. The efficiency of up-conversion of the 1 mm radiation into the third harmonic 
by a single-layer graphene is about 10 -7 at pumping level ~ 100 kW/cm 2 , room temperature, and 
concentration ~ 5 x 10 11 cm -2 . 

PACS numbers: 72.20.Ht, 72.80.Vp, 42.65. Ky 



The nonlinear response of carriers in graphene has been 
studied during the past years for the cases of the heat- 
ing by a strong dc electric field and of the interband 
pumping in the visible or IR spectral regions. In the mi- 
crowave frequency region, the nonlinear response is also 
effective because the massless and gapless energy spec- 
trum and a suppressed relaxation of low-energy carriers. 
[l| It is important that different semiconductor struc- 
tures is widely used for the multiplication (up-conversion) 
of the microwave pumping into THz radiation. [2| Re- 
cently the frequency multiplier in the GHz region, which 
is based on the nonlinear response of a bipolar graphene 
near the electroneutrality point, have been demonstrated. 
[H Such a multiplier should be effective in the frequency 
range u> <C v gr where v gr is the generation-recombination 
rate (in short devices, a fast contact injection can reduce 
this frequency restriction). For the high-frequency re- 
gion, the collisionless regime of the high-order harmon- 
ics generation was considered in Ref. 4. These calcu- 
lations are valid under the condition u S> v m which is 
correspondent to the THz pumping because the momen- 
tum relaxation rate v m ~ 10 13 s _1 for a doped sample 
with the mean momentum p. [l| The regime of nonlinear 
response under a strong pumping in the mm frequency 
range, when v m ^> cj, is not considered up to now. 

In this paper we examine the heating of carriers and 
the high-order harmonics generation in doped graphene 
under a strong ac electric field E t = E cos u>t of the mm 
frequency range, so that ui <C v m . We restrict ourselves 
by the weak anisotropy case, eE jv m <C p, when the elec- 
tron distribution function f pt + Af pt is written through 
the isotropic part, f p t, and the frequency- independent 
anisotropic addition 
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where v v is the momentum relaxation rate. The symmet- 
ric part of distribution is governed by the averaged over 



p-plane equation (see similar calculations for heating by 
a dc field in Ref. 5) 
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Here the collision integrals J r (ft\p) describe the nonelas- 
tic relaxation by acoustic phonons (r = ac) and the 
carrier-carrier scattering (r = cc). We neglect the re- 
laxation by optical phonons because the energy of hot 
carriers is less than the optical phonon energy. The Joule 
heating contribution is described by 



dAf pt _ (eE) 2 8 



dp 



2p dp 



cos 2 wt, (3) 



where a time dependency of f p t appears due to the pro- 
portional to cos 2 uit factor. 

The concentration of carriers and the current density, 
n and I t , are determined through f pt and Af pt according 
to the standard formula 



(4) 



where v p = vp/p is the velocity of electron with the 
momentum p and v ~ 10 8 cm/s is the characteristic 
velocity of the linear dispersion law. Using Eq. (1) and 
introducing the nonlinear conductivity at according to 
It = (TiEt, one obtains 
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The right-hand part here is transformed for the case of 
the short-range scattering, when v p ~ vap/h is written 
through the characteristic velocity v^, this approach is 
valid up to energies ~100 meV. p Here ft = f P =ot stands 
for the time-dependent maximal distribution at p = 0. 

Because of predominance of the intercarrier scatter- 
ing in a heavily doped graphene if w < v cc , the colli- 
sion integral J cc imposes the quasiequilibrium distribu- 
tion f pt = ft/[exp(vp/T t )(l-ft)+ft] written through the 
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time-dependent effective temperature, T t , and ft- These 
parameters are related by the normalization condition 
and the energy balance equation. Introducing the inte- 
grals over the dimensionless momentum 
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one transforms the normalization requirement in Eq. (4) 
as follows 
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The balance equation for the energy density, 
4 J dpf pt vp/(2Tr'h) 2 , is transformed into 
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FIG. 1: Time-averaged temperature T (a) and maximal dis- 
tribution / (b) versus pump intensity S plotted for concen- 
trations 5 x 10 11 cm" 2 (1) and 10 12 cm" 2 (2). 



d_ 

m 



hi > 

2 



fMf t ) 



v R 2 (T t y T? T t 



V TT \hvj h \ Iph J 

where the Joule contribution was written through the 
time-dependent conductivity a t oc /(, see Eq. (5). The 
energy relaxation contribution in the right-hand side of 
(8) is written for the case of the quasielastic scattering 
by acoustic phonons at temperature T p h- 0] Here we 
used the energy relaxation rate v e p/h written through 
the characteristic velocity v e <C v,Vd- 

Numerical solution of Eq. (8) is performed under the 
relation (7) and the periodicity conditions Tt_|_ w / W = T t 

and ft+n/u — ft- The nonequilibrium distribution f pt is 
determined through T t = T + 5T t and f t = f+6f t . Here 
we separated the time-averaged parts, T and /, which 
describe hot electrons, and the oscillating contributions, 
STt and 5ft- The calculations are performed at T p h =300 
K for the concentrations n = 5 x 10 11 cm" 2 and 10 12 
cm" 2 , the intensities up to 200 kW/cm 2 , and the char- 
acteristic velocities Vd ~ 3 x 10 r cm/s and v e ~51.4 cm/s. 
The dependencies of time-averaged contributions T and 
/ on the pumping intensity S are shown in Figs, la and 
lb. At low intencities of pumping S <5 kW/cm 2 , T t and 
ft vary strongly in analogy with to the results on heating 
by a dc electric field. @ A slowness of these dependencies 
takes place at higer pumping levels, in the region S ~10 - 
100 kW/cm 2 . Since the relaxation rates increase with en- 
ergy, both T t and ft increase with concentration. Notice, 
that a response on a probe dc field (photoconductivity) 
is described by Eq. (5) with the time-averaged / so that 
Fig. lb determines the photoconductivity versus S. 

The time-dependent contributions ST t = T t — T t= o & n d 
5ft = ft—ft—o are shown in Figs. 2a and 2b for the region 
S =25 - 100 kW/cm 2 at concentration 5 x 10 11 cm" 2 . 




FIG. 2: Nonharmonic time-dependent contributions to tem- 
perature STt (a) and to maximal distribution 5ft (b) for the 
concentration 5 x 10 11 cm" 2 and the 1 mm pumping intensi- 
ties 25 kW/cm 2 (1), 50 kW/cm 2 (3), and 100 kW/cm 2 (3). 



In particular, 5ft determines the time-dependent part of 
at, i- e. the nonlinear current under a strong pumping. 
Similar dependencies for the concentration 10 12 cm -2 are 
not shown. 

Further, we discuss the efficiency of transformation 
of the millimeter pumping into THz signal for the 
case of normal propagation both the pumping wave 
and the fcth order response determined by the field 



E z exp(— ikuit). 
density of Eq. 
I* = Sfe Ifcw ex P( — ikut) with the nonzero odd harmon- 
ics (it is due to in-plane isotropy of the problem). The 



Substituting f pt into the 2D current 
(4) one obtains the Fourier expansion 
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fcth harmonic of in- plane THz field, E z , is governed by 
the wave equation: [7] 

d 2 E z fkuj\ 2 4tt, v , , 

-jf + 'l-) Ez + i-^kuI ku 5 N {z) = 0, (9) 

where e is the uniform dielectric permittivity and 8n(z) 
is the transverse form-factor of the nonlinear current cal- 
culated above. For the TV-layer graphene stack, this fac- 
tor is normalized as J d ^ 2 dz5N(z) — N with d — > 0. 
The solution of Eq. (9) at z ^ takes form E 2 = 
Efc w exp(iKkui\z\) with the wave vector Kku = y/eku)/c. 
The amplitude of fcth harmonic, Efe w = — ATIfc w (27r/\/ec) 
is obtained from the boundary condition at z — ^ 0. 
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FIG. 3: Efficiency of transformation into fcth harmonics, 
Skui/S, versus pumping intensity S for fc = 3 (a) and fc = 5 
(b) at doping levels 5 x 10 11 cm -2 (1) and 10 12 cm -2 (2). 

The efficiency of up-conversion of 1 mm pumping of in- 
tensity S into THz signal of intensity Shu is determined 
by the ratio Skcu/S = 4|E fcw | 2 /£' 2 . Substituting I ku from 
Eq. (4) we plot this efficiency for fc =3 and 5, see Figs. 3a 
and 3b, where third and fifth harmonics are correspon- 
dent to 0.9 THz and 1.5 THz signals, respectively. One 



obtains the output intensities Sz u ~1.1 - 8.9 mW/cm 2 
for the pumping range 50 - 150 kW/cm 2 at concentration 
5 x 10 11 cm -2 while S^u corresponds up to ~10 fiW/cm 2 
range. Since Sku> oc A^ 2 for the Nth layer structure of 
epitaxial graphene, Q the up-converted 0.9 THz signal 
appears to be about 0.1 W/cm 2 (or ~0.2 mW/cm 2 at 1.5 
THz) under the 100 kW/cm 2 pumping at the 5 x 10 11 
cm~ 2 doping level and N =5. 

Next we briefly discuss the assumptions used in our cal- 
culations. We describe the heating of carriers with the 
use of the quasiequilibrium distribution approach. This is 
valid for the heavily doped graphene which is suitable for 
an efficient up-conversion. Typically, v m ~ v cc in doped 
samples and, together with the low-frequency condition 
u) < v m , we arrive at the restriction v e < uj < u m , v cc 
which is correspondent to the 10 9 s _1 - 10 13 s -1 frequency 
region. In spite of f p= ot < 1 we neglected the interband 
generation-recombination processes because the density 
of states vanishes at low energies. Also, the interband op- 
tical transitions are not essential in the spectral region up 
to bu. Beside of this, we perform the simplified electro- 
dynamical estimates of the efficiency of up-conversion for 
the case of the normal propagation of radiation through 
the structure with the same dielectric permittivities of 
the substrate and cover layers. The simplifications listed 
do not change either the character of the high-order non- 
linear response or the numerical estimates for the effi- 
ciency of up-conversion. 

In summary, we have demonstrated that the heating 
by a strong ac field results in an essentially nonhar- 
monical response. This leads to the effective generation 
of the third harmonic and opens a possibility for the 
efficient up-conversion of mm pumping into THz signal. 
An experimental study of the process considered and the 
electrodynamical calculations of a graphene structure 
integrated to an antenna or into a THz resonator are 
very opportunely now. 



The author is grateful to E. I. Karp for insightful com- 
ments. 



[1] A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. 

Novoselov, and A. K. Geim, Rev. Mod. Phys. 81, 109 

(2009); N. M. R. Peres, Rev. Mod. Phys. 82, 2673 (2010). 
[2] P. H. Siegel, IEEE Transactions on microwave theory and 

techniques 50, 910 (2002); A. V. Raisanen, Proc. IEEE 

80, 1842 (1992). 
[3] H. Wang, D. Nezich, J. Kong, and T. Palacios, Electron 

Device Letters, IEEE 30, 547 (2009). 
[4] S. A. Mikhailov and K. Ziegler, J. Phys.: Condens. Matter 

20, 384204 (2008). 
[5] O. G. Balev, F. T. Vasko, and V. Ryzhii, Phys. Rev. B 

79, 165432 (2009); O. G. Balev and F. T. Vasko, J. Appl. 

Phys. 107, 124312 (2010). 



[6] F. T. Vasko and V. Ryzhii, Phys. Rev. B 77, 195433 
(2007). 

[7] P. N. Butcher and D. Cotter, The Elements of Nonlinear 
Optics (Cambridge Univ. Press, Cambridge, 1990); F. T. 
Vasko and A. V. Kuznetsov, Electron States and Optical 
Transitions in Semiconductor Hetero structures (Springer, 
New York, 1998). 

[8] W. A. de Heer, C. Berger, X. Wu, M. Sprinkle, Y. Hu, 
M. Ruan, J. A. Stroscio, P. N. First, R. Haddon, B. Piot, 
C. Faugeras, M. Potemski, J.-S. Moon, J. Phys. D: Appl. 
Phys. 43 374007 (2010). 



